Emulating solid-state physics with a hybrid system of ultracold ions and atoms 
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We propose and theoretically investigate a hybrid system composed of a crystal of trapped ions 
coupled to a cloud of ultracold fermions. The ions form a periodic lattice and induce a band struc- 
ture in the atoms. This system combines the advantages of scalability and tunability of ultracold 
atomic systems with the high fidelity operations and detection offered by trapped ion systems, ft 
also features close analogies to natural solid-state systems, as the atomic degrees of freedom couple 
to phonons of the ion lattice, thereby emulating a solid-state system. Starting from the micro- 
scopic many-body Hamiltonian, we derive the low energy Hamiltonian including the atomic band 
structure and give an expression for the atom-phonon coupling. We discuss possible experimental 
implementations such as a Peierls-like transition into a period-doubled dimerized state. 
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The study of quantum-many body systems through 
quantum simulators, as initially proposed by Feynman 
[iHS], is one of the major challenges of modern ex- 
perimental quantum physics. Amongst others, trapped 
atomic gases cooled to quantum degeneracy and ions in 
electromagnetic traps have proven to be two of the most 
successful systems for this task [H [5]. These two sys- 
tems feature complementary strengths and advantages. 
Trapped ions allow for remarkably precise control and 
read-out of the vibrational and internal degrees of free- 
dom and can be used to construct a universal quantum 
simulator [6j. Their scalability to large systems however 
remains an important issue and is currently under intense 
investigation |7l|8]. Atomic systems, on the other hand, 
can easily incorporate large particle numbers. However, 
quantum operations can typically not be performed with 
very high fidelity in this case. Both systems are tun- 
able throughout a wide range of parameters and have 
recently been used in a number of seminal quantum sim- 
ulations [HIS]. 

In this Letter we propose a setup consisting of the 
combination of a trapped ion crystal interacting with a 
system of ultracold fermionic atoms, combining the ad- 
vantages of both systems [9-15 . We envision fermionic 
atoms simulating electrons in a solid-state system with 
ions forming the lattice. A crucial component of solid- 
state systems, namely the phonon-electron coupling, 
arises naturally in our scenario. This feature is absent in 
ultracold atoms trapped in optical lattices, where there 
is no significant back action of the atoms on the lattice. 
Likewise, trapped ion systems do not naturally feature 
fermionic statistics, so that the quantum simulation of 
coupled fermion-phonon systems would require a digi- 
tal decomposition into many quantum gates [16j adding 
significant experimental complexity. Besides the atom- 
phonon coupling and fermionic statistics, our system in- 
cludes a number of genuine solid-state effects from the 



outset such as sharply peaked lattice potentials, finite- 
size effects and spin dependence. The proposed simu- 
lator may be used to study models of electron-phonon 
coupling, such as the Frohlich model which describes oc- 
currence of polaronic quasiparticles, Peierls transitions 
and the emergence of phonon-mediated interactions. 

Alternative proposals for investigating phonon cou- 
plings include self-assembled lattices of ultra-cold dipolar 
molecules |17l and optically trapped nano particles '18], 
whereas back-coupling to cavity photons has also been 
demonstrated |19j . Our proposal is complementary to 
these approaches and offers significant advantages in the 
ability to control precisely the couplings and configura- 
tion of the system. 

The interaction between an atom and an ion is caused 
by an induced dipole moment in the atom, and for 
sufficiently large distances is given by VAi{r) — ^ 
where r denotes the distance between atom and ion 
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Solid state ''Li-i^^Yb+ *°K-«Ca^ 



Lattice spacing d (nm) 

Length scale R* (nm) 

Energy scale E* (kHz) 

d/R* 

mi/ruf 

Fermi energy (MHz) 

Phonon energy (MHz) 



0.3-0.6 

0.026 

10" 

10-20 

lO^-lO^ 

10* 

10'' 



lO'^-lO* 

71 

166 

14-140 

29 

0.02 

0.01-10 



lO^'-lO* 

245 

2.1 

4-40 

1.0 

0.02 

0.01-10 



TABLE L Comparison between a typical natural solid-state 
system and the proposed quantum simulator. The charac- 
teristic length and energy scales for atom-ion interaction are 
defined in the text and correspond to the hydrogen electron 
for the solid-state system. The mass of the fermion (elec- 
tron/atom) and ionic core are denoted by m/ and mt respec- 
tively. A mono-atomic crystal of elementary atoms (e.g. Na) 
is assumed for the solid-state system. 



ionic phonon dispersion 



atomic dispersion 
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FIG. 1. A cloud of ultracold atoms is brought close to a linear crystal of ions, which are trapped in a Paul trap (a). Due to 
the strong Coulomb repulsion, the ions are strongly pinned to a periodic lattice and the excitations are bosonic phonons, with 
a typical spectrum ranging over a few MHz shown in (b). To lowest order and for transversal atom-ion separation Ar — 0, 
mobile atoms experience a periodic spatial potential via the attractive atom-ion potential, forming Bloch states (shown in (d)) 
and energy bands (c) for d = 15 R* and ID scattering phases ipo = —ifie = 7r/4. 



This interaction defines both a natural length scale 
i?* = ^J^hriAC^J^ ^ which is typically in the 100 nm 
regime [20], and an energy scale E* = /i*^/(4m^C4), 
where m^ is the atomic mass [3T]. The ions repel each 
other with a long range r~^ Coulomb interaction leading 
to phonon branches in the ion string which in turn couple 
to the atoms. Interactions between the atoms arc short- 
range and can be tuned e.g. by Feshbach resonances. In 
Table 1, we list some typical energy and length scales of 
our proposed quantum simulator as compared to typical 
natural solid-state systems. 

One advantage of our proposal is that the mass ratio 
between the lattice species and the fermionic species can 
be chosen to emulate different phenomena. The large 
mass ratio of ^Li atoms and ^^''Yb"'" ions is reminiscent 
of solid-state systems, whereas ^°K and ''*'Ca+ involves 
nearly equal masses. Alternatively, the combination of 
^°K and ^Be+ allows for the consideration of an inverted 
mass ratio regime and emphasizes the motion of the ions. 

The outline of this Letter is as follows: we first derive 
the effective low energy Hamiltonian of the system of ions 
and atoms, both occupying a single spin state, and de- 
termine an explicit expression for the atom-phonon cou- 
pling. Thereafter we demonstrate the use of our sys- 
tem as a quantum simulator by investigating a Peierls- 
like transition of a system of fermionic atoms. Finally, 
we discuss the experimental implementation and further 
possibilities in an outlook. 

Derivation of the low energy model The system is mi- 
croscopically described by three terms in the Hamiltonian 
T-L =%! + Ha + 'Hai, where 
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describes A^i ions, each of mass rrii and charge -\-e trapped 
in a Paul trap represented by the single-particle poten- 
tial V/(R„), with pairwise Coulomb interactions. R„ 
and P„ are the position and momentum operators of 
the n-th ion respectively. The A^a atoms are described 



by a Hamiltonian "Ha = J2 
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given by VA(r) which is oriented parallel to, and can be 
displaced by a distance A^ from, the center of the Paul 
trap. In contrast to the case of ions [22], the fermionic 
quantum statistics of the atoms plays an essential role. 
The atom-ion interaction is described by the Hamiltonian 
Hai = J2i j ^Ai{\^i — RjD- For the theoretical descrip- 
tion, we consider an infinitely long system or equivalently 
a system with periodic boundary conditions. To derive 
a low energy effective Hamiltonian, we separately trans- 
form the ionic and atomic degrees of freedom into diag- 
onal form and express the coupling between the systems 
in this basis. 

The dominant energy scale is set by the interionic re- 
pulsion, which, in conjunction with the trapping poten- 
tial Vi and a given axial ionic lattice spacing d (typically 
of the order l-20fim) defines the classical equilibrium po- 
sition of the j-th ion Rj. This classical configuration 

minimizing the electrostatic energy Ej corresponds to 
a linear string, planar zigzag or structures of three di- 
mensional crystals \23, . As the fluctuations of quantum 
or thermal origin in the ion displacement are small, the 
ionic Hamiltonian is well suited to an expansion in terms 



of the position fluctuations operators STlj 



Hence we write Hi = E 
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Rj Rj. 

. . . where 



the first order terms vanish identically and the low en- 
ergy spectrum and the dynamics are determined by the 

structure of ^f - E„ A?/2m/ + | E„,„' Vn,n'SRnSRn', 
where we use the combined label n = (j, D) for the 
coordinate of the j-th ion in spatial dimension D and 

Vr,,n'=d''Vi/d5Rr,d6Rr,,. 

After transforming into the basis of phonons with 
bosonic operators [as,Q;J] — 6s^s', where s is a collective 
mode label |24|, the Hamiltonian takes on the form 



Hi = i?}°) +^^,(ata, + -) + — +0(«3), (2) 



2m 



eludes a tight radial optical dipole (or magnetic) trap 



where V is an effective momentum operator associated 
with a free particle of mass rh (see Appendix [b]) and we 
have set H = 1. The spectrum ujg for a particular linear 
ion string is shown in Fig. [1]d. 

In a similar fashion, an expansion in SRj (for para- 
metrically fixed r^) can also be made for the atom- ion 



operator Hai — T~i 
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To lowest order, "Haj, the atoms feel a periodic lattice 
of ions pinned to their classical equidistant positions and 
the higher orders CHajj etc.) describe couphngs to ionic 

phonon modes. The atomic eigenstates of T-La +'^ai ^^^ 
Bloch states 0fc_Q(r) characterized by a quasi-momentum 
k and a band index a. The form of the long-range atom- 
ion potential VAi{r) is unbounded from below and the 
low-lying energy bands (see Fig.nt ) are essentially super- 
positions of local bound states of each ion with flat dis- 
persions (bandwidth < 20 Hz for Yb+-Li with d = 1 /im) . 
Furthermore, they are energetically separated by a large 
gap (see Fig. flfc) and should not play a significant role if 
the atoms are moved into the ion string in an appropriate 
sequence [25]. On the other hand, the Bloch states with 
energies above the threshold (see upper state in Fig. [l];) 
have a considerably larger bandwidth of typically a few 
kHz. We expect that adiabatic loading of the atoms into 
the ion system will populate the lowest of these bands. 
As it is well separated in energy from other bands, we can 
work within a single-band approximation and henceforth 
omit the band index. Alternatively, the system may be 
prepared into this band by either cooling the atoms or 
the ions as long as the coupling to the molecular bands 
is low. When the atom and ion traps are separated by a 
large distance A^ the atomic band structure is bounded 
from below and it is most relevant to consider the band 
with lowest energy. 

The single-atom part of the Hamiltonian is 
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where, working within the atomic lowest band approx- 
imation, cj. — J d'^r <l)k{r)'^^{r) creates an atom in the 
Bloch mode k. As the long-range atom- ion coupling di- 
verges for small separations, an effective description of 
the short-range coupling is important. Although an exact 
calculation of the 3D atomic Bloch states cannot be per- 
formed by typical band structure calculations, which gen- 
erally require Bloch waves to be defined over all space, it 
is possible to determine the effective ID Bloch states and 
the corresponding dispersion relation within quantum de- 
fect theory, as discussed in Appendix |A] For sufficiently 
strong radial confinement, any low energy eigenstate will 
be well approximated by the form ipiz, p) = tpioiz) 4'i-{p) 
in spatial regions between the ions, where the potential 
felt by the atoms is almost cylindrically symmetric. Here 
4'±{p) is the lowest 2D harmonic oscillator state of the 
optical dipole trap and the trap is oriented in the e^- 
direction. The effect of the ion on the wave function can 
then be characterized by phase shifts given by the odd 
and even ID phases ipo and <pe |20j . A typical dispersion 
relation and Bloch state and are shown in Fig. [l]; and 



[TJi respectively. The effective ID theory is valid only for 
atom-ion distances that are larger than the radial wave 
packet size l±. For Li and Yb+, and assuming a radial 
atomic trap frequency of ioj_ — 2'k x 200 kHz, we find 
l± ^ R* . Hence for these parameters, the ID description 
is applicable only for inter-ion distances d ^ R* and en- 
ergies E <^ E* . For tighter atomic confinement or heav- 
ier atomic species, the description holds over larger pa- 
rameter regimes. Obtaining the exact form of the atomic 
wave function closer to the ion would, however, require a 
3D description. 

The dominant term of the coupling Hamiltonian is 
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which preserves the atom number but couples states of 
differing phonon number. It is useful to define the po- 
tential 
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(r - Rj) ■ ep 
Ir-R,|6 ■ 
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where i — i{n) and D = D{n) refer to the ion number and 
spatial dimension. In terms of this potential the phonon 
coupling is given by 
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where the quantities u„ , Vn and r2„ are determined in 
the phonon transformation (see Appendix IbI) . Note that 

^k k' ^^ only uniquely defined up to the (arbitrary) com- 
plex phases of the individual Bloch and phonon modes. 
To calculate A^ ^, in the case of overlapping atom and 
ion traps, we require an accurate description of the Bloch 
states very close to the ions, which is unfortunately not 
available in the quantum defect calculation. Further- 
more, the secular approximation for the ions neglecting 
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FIG. 2. The interaction between atoms and ions is expressed 
as an atom-phonon coupling \^l^^, '^'^'| in the low energy 
limit. For atoms displaced by Ar = 0.75 R* and tuning to 
a small oj^^^/tj .^ = O.OOQSJli; = 30 -E*, the longitudinal cou- 
pling (a) is several orders of magnitude smaller than the trans- 
verse coupling (b) and the former may exceed the band width, 
while still being smaller than the band gap. 
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FIG. 3. Experimentally tuning the ion trapping allows for the observation of a Peierls-like phase transition when atoms are 
confined in an effective ID dipole trap at a separation Ar from the ion string. The discrete translational symmetry of the linear 
chain (a), is spontaneously broken towards a zigzag pattern of the ions (c) and a band gap opens at fc = 7r/2d for the atoms. 
At different atomic fillings n (per ion), and below a critical temperature Tc, a phase transition occurs (b). Unlike the Peierls 
effect in solid-state materials, the optimal filling that generates the largest Tc can be larger than n = 1/2 (d). 



ion micromotion may break down [26L I27| . We note that 
a full quantum mechanical treatment of the effect of the 
time-dependent trapping potential will require a many- 
body quantum Floquct formalism [2 8) and is beyond the 
scope of the present paper. 

For an atomic cloud tightly confined in the radial direc- 
tion and radially displaced from the ion string by Are^, 
the low energy single particle Bloch states are in the ra- 
dial ground state and of the form ipk^{r) = 'ipk^iz)(j>±{p)- 
The ID Bloch (f>k{z) states in the potential Vi]j{z) — 
X^i VAiilzez — Hil) can be calculated by diagonalizing the 
Hamiltonian in the plane wave basis. The atom-phonon 
coupling A^**^, , which vanishes unless q ^ k — k' where q 
is the quasi-momentum corresponding to phonon index 
s, has been calculated (see Appendix [C|) and is shown in 
Fig. [2J The transversal coupling features a pronounced 
peak at g = |fc — A:'| = Tr/2d as the system approaches the 
zigzag transition, discussed in the following. 

Transverse Peierls-like transition As a demonstra- 
tion of our quantum simulator, we investigate a Peierls- 
type transition [21]. In solid-state a metallic ID crys- 
tal can undergo a transition to an insulating phase, the 
so-called Peierls transition, when cooled below a critical 
temperature. This is because an increased period of the 
ion lattice is energetically favorable since this allows a 
gap to open in the electronic band structure. For par- 
ticular fermionic band filling, this can decrease the elec- 
tronic energy enough to offset the energy increase of the 
ion distortion. In our case, the energy that the linear 
ion string gains as it is distorted into a period-doubled 
zigzag configuration is offset by the lowering of the en- 
ergy of the atomic system at half- filling (see Fig. Isfc). 
Analogous distortions can be observed at other rational 
fillings. 

The structural phase transition from a linear to a 
zigzag configuration in trapped ions has previously been 
the subject of theoretical study as it can be coupled to 
quantized fields [30H32] . The present proposal goes one 
step further by using a fermionic quantum field to induce 
the phase transition, thereby emulating the Peierls effect 



as observed in natural solid-state systems. 

We investigate the effect with a separation A^ be- 
tween the atom and ion traps and treat the ion posi- 
tions classically. We perturb the ions from their equi- 
librium positions by an amount S oriented along the 
trap separation such that R^ = R^ -t- (—1)"^. The en- 
ergy increase of the ions due to the zigzag configura- 
tion AEi = miu;l^^^^jy2 + 0{S*), with S = \5\, is 
associated to a zigzag oscillation mode with frequency 
L^k=iT/d,x- For T — >■ 0, zigzag order will always be pre- 
ferred as the band gap scales linearly with S and the 
energy of classical ion displacement increases with S^ , for 
small S. At finite temperature, the preferred configura- 
tion is determined by minimizing the free energy. To this 
end, we calculate the free energy of the non-interacting 
atomic system by parametrically varying 6 and subtract 
the energy gained by the ions in the presence of the trap. 

A sample band structure for non-interactingj"'K atoms 
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near a ^"^Ca"*" ion string is shown in Fig. pp where a 



zigzag displacement S 



u = {n. 
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(0) 



0.019 R* minimizes the energy 
) of a system with atomic half- 



filling and static ion positions Rj with uJzz ~ 27rx 60 kHz, 
ion spacing d = 2 fim w 8 i?* and atomic temperature 
Ta = 15 nK w 0.15 E*/kB. The phase diagram for 
non-interacting '*°K with '^^Ca^ is shown in Fig. [Sp for 
varying temperature and atom-ion trap separation. Fur- 
thermore, while the largest critical temperatures occur 
for half-filling when A^ is large, the optimal critical tem- 
perature, shown in Fig. [Sli, occurs for larger fillings per 
ion, n, as A,, — > 0. 

Implementation For the proposed quantum simula- 
tion, the interacting quantum systems can be initialized 
individually and their mutual coupling strength can be 
controlled by varying the distance A^ between ions and 
fermions. The ion crystal can be cooled to its ground 
state using resolved sideband cooling and in order to 
reach the regime in which the low frequency phonon 
modes can be resolved, one may cool in a tight trap and 
then adiabatically reduce the trap frequency |33) . The 
motional state of the ion crystal can be initialized and 



read out by laser-induced coupling to internal states of 
the ions which can be measured with standard fluores- 
cence measurements [51]. Building on these techniques, 
interferometric methods have been proposed to probe ion 
crystals near structural phase transitions [55!. 

The possibility of controlling small numbers of ultra- 
cold fermions in optical traps has recently been demon- 
strated |5nj- The atomic density profile is accessible via 
absorption imaging in combination with time-of-flight 
analysis. For low densities, we can retrap the atoms 
in an optical molasses trap to prevent them from leav- 
ing the imaging region during detection. While we have 
focused on fermionic atomic species in this paper, it is 
also possible to trap bosonic species such as ''Li, where 
ion-controlled tunneling dynamics and Josephson physics 
may be investigated |57| . 

Hybrid ultracold atom-ion systems pose several chal- 
lenges for experimental realization. Inelastic collisions - 
in particular charge transfer and molecular ion formation 
- may limit the lifetime. It has been observed that for 
some species the rates may be small enough to study in- 
tricate setups, such as in a recent experiment where the 
spin relaxation of a single Yb+ ion immersed in a bosonic 
Rb bath has been studied [T5]. However, no experimental 
data are available for inelastic rates between ultracold K 



or Li with trapped ions. The use of polarized fermionic 
systems may reduce the effects of three-body recombi- 
nation as binary atom collisions are suppressed in the 
ultra-cold regime [38 . In the geometry in which we stud- 
ied the Peierls transition, ions and atoms are separated 
and inelastic collisions are suppressed. For the reduction 
of atomic heating effects due to ion micromotion we rely 
on the largest experimentally feasible mass-ratio between 
ions and atoms (^^^Yb"*" and ^Li) [57j and preventing the 
atoms to approach the ions too closely [26]. To further 
mitigate micromotion-induced heating, it may be feasible 
to use even heavier molecular ions, which could be sym- 
pathetically cooled by suitable co-trapped ion species. 

Future work will focus on solving the 3D many- 
body atom-ion problem including short-range interac- 
tions, considering the time-dependent potential of the 
ions and the role of the spin degrees of freedom in order 
to engineer entanglement between the atoms and ions. 

This work was supported by the DFG via Son- 
derforschungsbereiche SFB-TR/49 and SFB-TR/21, 
Forschergruppe FOR 801, the excellence cluster "The 
Hamburg Centre for Ultrafast Imaging" of the DFG 
and the excellence cluster under contracts PICC and 
AQUTE. RG acknowledges fruitful discussions with 
S. T. Dawkins and K. Jachymski. 



[1] R. Feynman, Int. J. Theoret. Phys. 21, 467 (1982). 

[2] S. Lloyd, Science 273, 1073 (1996). 

[3] J. I. Cirac and P. Zoller, Nat. Phys. 8, 264 (2012). 

[4] I. Bloch, J. Dalibard, and S. Nascimbene, Nat. Phys. 8, 

267 (2012). 
[5] R. Blatt and C. F. Roos, Nat. Phys. 8, 277 (2012). 
[6] B. Lanyon ei al, Science 334, 57 (2011). 
[7] R. Islam ei al, arXiv:1210.0142 (2012). 
[8] J. Britton et al, Nature 484, 489 (2012). 
[9] A. T. Grier, M. Cetina, F. Orucevic, and V. Vuletic, 
Phys. Rev. Lett. 102, 223201 (2009). 
[10] C. Zipkes, S. Paltzer, C. Sias, and M. Kohl, Nature 464, 

388 (2010). 
[11] S. Schmid, A. Harter, and J. H. Denschlag, 

Phys. Rev. Lett. 105, 133202 (2010). 
[12] C. Zipkes ei al, Phys. Rev. Lett. 105, 133201 (2010). 
[13] W. G. Rellergert ei al, Phys. Rev. Lett. 107, 243201 

(2011). 
[14] L. Ratschbacher, C. Zipkes, C. Sias, and M. Kohl, Nat. 

Phys. 8, 649652 (2012). 
[15] L. Ratschbacher ei al, arXiv: 1301. 5452 (2013). 
[16] J. Casanova, A. Mezzacapo, L. Lamata, and E. Solano, 

Phys. Rev. Lett. 108, 190502 (2012). 
[17] G. Pupillo ei al., Phys. Rev. Lett. 100, 050402 (2008). 
[18] M. GuUans ei a/., Phys. Rev. Lett. 109, 235309 (2012). 
[19] R. Mottl ei al., Science 336, 1570 (2012). 
[20] Z. Idziaszek, T. Calarco, and P. Zoller, Phys. Rev. A 76, 

033409 (2007). 
[21] Note that the common definition of R* and E* makes use 
of the reduced mass of the two-body system of an atom 
and an ion. However, our definition is better suited to a 
many-body system. 



[22] Due to the very strong interactions and large interparticle 

spacing, the ions have a very low phase space density. 
[23] D. F. V. James, Appl. Phys. B 66, 181 (1998). 
[24] For instance, in the linear crystal regime s = {k, 7) and 

7 denotes either the acoustic or one of the two transverse 

optical modes. 
[25] H. Doerk, Z. Idziaszek, and T. Calarco, Phys. Rev. A 81, 

012708 (2010). 
[26] L. H. Nguyen, A. Kalev, M. Barrett, and B.-G. Englert, 

Phys. Rev. A 85, 052718 (2012). 
[27] M. Cetina, A. T. Grier, and V. Vuletic, Phys. Rev. Lett. 

109, 253201 (2012). 
[28] H. Landa, M. Drewsen, B. Reznik, and A. Retzker, New 

J. Phys. 14, 093023 (2012). 
[29] R. Peierls, Princeton Series in Physics (1991). 
[30] E. Shimshoni, G. Morigi, and S. Fishman, 

Phys. Rev. Lett. 106, 010401 (2011). 
[31] C. Cormick and G. Morigi, Phys. Rev. Lett. 109, 053003 

(2012). 
[32] J. D. Baltrusch ei al, Phys. Rev. A 84, 063821 (2011). 
[33] G. Poulsen and M. Drewsen, arXiv: 1210.4309 (2012). 
[34] D. Leibfried ei al, Phys. Rev. Lett. 77, 4281 (1996). 
[35] G. D. Chiara, T. Calarco, S. Fishman, and G. Morigi, 

Phys. Rev. A 78, 043414 (2008). 
[36] F. Serwane ei al.. Science 332, 336338 (2011). 
[37] R. Gerritsma ei al, Phys. Rev. Lett. 106, 060503 (2011). 
[38] A. Harter ei al, Phys. Rev. Lett. 109, 123201 (2012). 



Appendix A: ID band structure within quantum 
defect theory 



Here, the short range behavior of the atom- ion interac- 
tion is described by quantum defect parameters and we 
work in the ID approximation [5D] . For this to apply, all 
energy scales in the system have to be small compared to 
the radial oscillator frequencies, such that the probability 
for any atom to be in a radially excited state is strongly 
suppressed. In this situation the Schrodinger equation in 
units of R* and E* for a single atom interacting with an 
ion located at R = for z ^ l± = ^/i/(2m^w^^j,) is 
given by 
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Here we assumed that the ionic mass is much larger than 
the atomic mass, mj >> m^, such that m^ « /i, the 
reduced mass. For a given finite energy E, this can be 
neglected in the limit of small atom-ion separation z — >■ 0, 
and the two-dimensional solution space approaches the 
asymptotic solutions 
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which are characterized by the even and odd phases (pe 
and (fio- These effective ID phase shifts fully character- 
ize the short-range interaction within the ID regime and 
have a non-trivial dependence on the 3D s-wave scat- 
tering length. Within quantum defect theory, an en- 
ergy eigenfunction is described by a superposition of the 
asymptotic solutions Eq. ( A2 ) and ( A3 ) in the vicinity 
of each ion up to a sufficiently small distance r^ 3> l±. 
Outside Tg, the two solutions ipeiz) and ^po{z) are propa- 
gated numerically up to the border of the real space unit 
cell for a fixed set of {E,ipe,^o) with the initial condi- 
tions chosen to match the asymptotic solutions smoothly 
at re- The two symmetric and anti-symmetric solutions 
span the corresponding solution space, guaranteeing that 
the Bloch state can be written as 
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with coefficients cl and cj, for the even and odd states 
respectively. Furthermore, in accordance with the Bloch 
theorem, each Bloch state fulfills 
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under spatial translation by a lattice distance d. If E 
does not lie within a band gap, the associated quasi- 
momentum k is defined by the constraint of both the 
wave function and its derivative being continuous at the 
edge of the unit cell. If the ion is chosen to lie in the 
middle of the unit cell at z = 0, symmetry implies 



M-l) = Mi): M-i) = -Ml), M-i) = -Mi), 

^q(— |) = i^'oii) and the matching conditions imply the 
linear relations 
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where 
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For normalizable, non-trivial solutions of (ce tCo ) to 
exist, the determinant of the coefficient matrix A(k) has 
to vanish. Thus, the quasi-momentum corresponding to 
the chosen energy E is efficiently determined using a root 
finding algorithm on the function (lei{A{k)). The Bloch 
states on the entire lattice are subsequently given by com- 
bining them piecewise in every real-space unit cell. 

For a single atom-ion system, the energy eigenstates 
can be classified into a set of discrete bound states (de- 
caying exponentially in real-space for large distances) and 
a continuum of unbound scattering states. In the limit 
of large intra-ionic spacings d ^ R* the Bloch states 
converge to these bound and scattering states and the 
band index takes on the role of the local eigenstate label. 
Only the bands corresponding to scattering states would 
have a significant band width (or correspondingly, sig- 
nificant hopping elements within the Wannier represen- 
tation) and the Bloch functions are exponentially sup- 
pressed in the regions between the ions. With decreasing 
d (or fixed d and decreasing atomic mass tua) the clas- 
sification into bound and scattering states loses its ap- 
plicability and also the Bloch states emerging from the 
bound states obtain a finite band width. 



Appendix B: The ionic phonon transformation 

Here we discuss the generalized unitary transforma- 
tion into the basis of phonon modes in both the linear 
ion chain, as well as the zigzag regime. In the low en- 
ergy regime the higher order terms are small and the low 
energy excitations are determined by the collective mode 

(2) 

structure of nj . Instead of obtaining the eigenmodes in 
the classical limit, it is favorable to perform the transfor- 
mation to the phonon basis on an operator level, which 
yields the precise transformation relations and avoids the 
problems of ordering of higher order terms before quan- 
tizing the bosonic fields. Recalling that the matrix ele- 
ments are Ki,n' = d'^Vi/dSRndSRn', and defining the lo- 
cal harmonic oscillator frequencies J7„ = \/Vn,nlfni and 
the associated annihilation operator for each ion and spa- 
tial dimension &„ = •\/7nir2„/2[^A„ -\- iPn j (mSln)\ (here 



h = 1), the ion Hamiltonian can be written as 
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Here, b (b^') denote the column vectors of aU annihi- 
lation (creation) operators respectively and the phonon 

coefRcient matrix is of the form i/pi, = a * i * I with 

the 3A^/-dimensional submatrices A„^„' = 2rn'n "'"' 

and hn^n' = ^n,n' ^n + A„^„'. The phonon spectrum and 
operators are obtained by a symplectic diagonalization, 
which preserves the bosonic commutation relations of the 
phonon operators and reduces to a canonical transforma- 

_y(s) ) IS an eigen- 
vector obtained from the diagonalization to the positive 
phonon energy w^, the associated phonon annihilation 
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tion in the classical limit. If x'"^ 



operator is defined by as 
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vectors, where E 



defines the metric for this 



Un and Vn , appearing in Eq. (l7|, are then the n-th ele- 
ments of the subvectors u^**) and v^^^ respectively. How- 
ever, only 3Nj — 1 phonon modes exist, i.e. the matrix 
ETJph only possesses 3-/V/ — 1 linearly independent eigen- 
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symplectic (bosonic) space. The missing eigenvector is 
fundamentally related to the existence of a Goldstone 
mode and can only occur if there is one zero eigenvalue. 
Let p be the eigenvector of -ffph to an eigenvalue zero. 
Its norm vanishes with respect to the norm defined by E 
and we normalize it with respect to the euclidean norm 
p^^p = 1. To span the entire subspace, we additionally 
define the vector q, which also has a vanishing E-norm 
q^'Sq = 0, is orthogonal to all eigenvectors x*^") and y^*^ 
with respect to E and furthermore fulfills q^Ep = i. The 
scalar products 
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defines a generalized collective momentum and position 
operator respectively. Furthermore the relation qtEp — 
i- defines a generalized mass m, appearing in the Hamil- 
tonian in Eq. [5] Since the eigenvectors and q span the 
entire space, the phonon operators {a^}, {at} together 
with V and Q span the same operator space as the orig- 
inal {6„} and {b^^}. Thus any operator acting on the 
Hilbert space of the ions can be expressed exactly in 
terms of the former and in this sense Eq. [2] is the ex- 
act expression of the full many-body ion Hamiltonian to 
second order, requiring the term ■p^/(2m). 

We now explicitly consider the phonon structure and 
spectrum for the linear ion chain in the limit of a large 
system. Here, the second order Hamiltonian can be ex- 



pressed in terms of the creation and annihilation oper- 
ators ak,o = N^ ' ^ e~'''"^^bi^o and their conjugates 
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The 2x2 coupling matrices for the longitudinal and 
transversal case are given by 



jjk,x/y 
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respectively. Here we defined 
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where ^(s) is the Riemann zeta function, Lis(a::) is the 
polylogarithm. 

The eigenmode energies are found to be 
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for the longitudinal and transverse phonon branches re- 
spectively, where w/ is the transverse ion trapping fre- 
quency. If the two-dimensional eigenvector of the matrix 

jjk,n jg denoted by ( ' I , the explicit form of the 
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phonon operators is given by 



a^k,n = uljj a^k,D + Wfe,_D 4,D 



For an isolated ion system, the instability towards the 
formation of zigzag pattern (i.e. a macroscopic occupa- 
tion of the 7r/d transverse phonon mode) appears if the 
expression in the root of Eq. (B7| becomes negative, i.e. 
approximately for 
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Appendix C: Calculation of the atom-phonon 
coupling for displaced atoms 

For the case of the atom cloud's radial confinement 
being much smaller than the displacement A^, the inte- 
grand is only significantly non-zero on the tube where the 
atoms are located (enforced by the radial density |(/>j.(r)p 
appearing as a factor) and the effective potential becomes 



a function of z only: Ki(r) = Vn{z). For the longitudinal 
phonon mode D(n) = z this is of the explicit form 
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Without loss of generality, we have assumed that the sep- 
arations between the traps is in the x-direction and for 
D[n) = X we have 
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and otherwise Va 
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0. Here we define the function 
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, z for the trans- 



verse and longitudinal modes respectively, as well as the 
periodic function FF{z) = e^^''^ Ff'{z), which fulfills 

F^{z + d) — F^{z). Expressing each Bloch state as 
(f>k{z) = e*'^^Ufe(2:), where Uk{z) is also a periodic func- 
tion Uk{z) = Uk{z + d), the phonon coupling becomes 
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where the function BZ{p) folds any momentum vector 
p back into the first Brillouin zone by adding an integer 
multiple of the reciprocal lattice vector 2Tr/d. The Kro- 
necker delta reflects the overall quasi-momentum conser- 
vation of the combined atom-ion system. Since the Bloch 
states are normalized with respect to integration over the 
entire lattice consisting of Nj sites, the second line of 
Eq. ( |C3[ ) rapidly converges to a value independent of Nj, 
as the lattice size is increased. Note that whereas a Bloch 
state (pkiz) is independent of the Brillouin zone k is in 
4>k{z) — 4>k+2Tr/d{z), this is uot the case for the function 

Uk{z). In principle, an additional term ^j. j., AJ. j,, Qcyf^ c^,, 
appears in Eq. ([5]) , which corresponds to a coupling of the 
atomic system to the center of mass motion of the ions. 
However, a calculation shows that the coefficient A^"^,, 
vanishes and hence we exclude it from the low energy 
Hamiltonian. 



